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Abstract 

We have studied the phase diagram for chromatin within the frame- 
work of the two-angle model. Rather than improving existing models 
with finer details our main focus of the work is getting mathematically 
rigorous results on the structure, especially on the excluded volume 
effects and the effects on the energy due to the long-range forces and 
their screening. Thus we present a phase diagram for the allowed 
conformations and the Coulomb energies. 
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1 Introduction 



Each protein aggregate together with its wrapped DNA comprises a nucleo- 
some core particle with a radius of about 5nm and a height of about 6nm; 
with its linker DNA it is the fundamental chromatin repeating unit. It carries 
a large electrostatic charge p. Whereas the structure of the core particle 
has been resolved up to high atomic resolution j2], there is still consider- 
able controversy about the nature of the higher-order structures to which 
they give rise. When stretched the chromatin string appears to look like 
"beads-on-a-string" in electron micrographs. 

The recent accumulation of evidence that the chromatin structure above 
the level of the core particle plays a key role in determining the transcrip- 
tional status of genes and genetic loci EJ illustrates the critical importance 
of understanding the fundamental folding properties of nucleosome arrays: 
Studies of chromatin compaction in response to changes in the ionic envi- 
ronment |3] have established that the phenomenon can be accounted for by 
electrostatic interactions between DNA, histone proteins and free ions jH]. 
Major contributions to these interactions are provided by the N-terminal do- 
mains of the core histones, which contain roughly half of the basic amino 
acids of the octamer, and the C-terminal domains of the linker histones, 
which contain about | of the positive charges in these molecules. Indeed, 
chromatin compaction requires the presence of the core histone N termini 
and 30nm chromatin fibers are not formed in the absence of linker histones. 
However, the precise interactions that lead to specific chromatin higher-order 
structures have remained elusive. 

The beads-on-a-string structure can be seen clearly when chromatin is 
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exposed to very low salt concentrations, and is known as the 10-nm-fiber, 
since the diameter of the core particle is about lOnm. With increasing salt 
concentration, i.e. heading towards physiological conditions (c ^lOOmM), 
this fiber appears to thicken, attaining a diameter of 30nm. The absence of 
the extra linker histones (HI or H5) leads to more open structures; so it is 
surmised that the linker histones act near the entry-exit point of the DNA; 
they carry an overall positive charge and bind the two strands together lead- 
ing to a stem formation [7J |HJ E] • Increasing the salt concentration decreases 
the entry-exit angle a of the stem as it reduces the electrostatic repulsion 
between two strands. 

Longstanding controversy JTUJ [TTJ [T3] surrounds the structure of the 
30nm fiber, for which there are mainly two competing classes of models: the 
solenoid models f3J [H| and the crossed-linker- models ^HJ IH3 E|- m the 
solenoid model it is assumed that the chain of nucleosomes forms a heli- 
cal structure with the axis of the core particles being perpendicular to the 
solenoidal axis (the axis of an octamer corresponds to the axis of the superhe- 
lical path of the DNA that wraps around it). The DNA entry-exit side faces 
inward towards the axis of the solenoid. The linker DNA is required to be 
bent in order to connect neighboring nucleosomes in the solenoid. The other 
class of models assumes straight linkers that connect nucleosomes located on 
opposite sides of the fiber. This results in a three-dimensional crossed- linker- 
pattern. Such an arrangement with peripherally arranged nucleosomes and 
internal linker DNA segments is fully consistent with observations in intact 
nuclei and also allows dramatic changes in compaction level to occur without 
a concomitant change in topology. 

Images obtained by electron cryomicroscopy (EC-M) should in princi- 
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pie be able to distinguish between the structural features proposed by the 
different models mentioned above. The micrographs show a crossed- linker- 
pattern at lower salt concentrations and they indicate that the chromatin 
fiber becomes more and more compact when the ionic strength is raised to- 
wards the physiological value [7|. However, for these denser fibers it is still 
not possible to detect the exact linker geometry. Experiments on dinucleo- 
somes have been made to check whether the nucleosomes collapse upon an 
increase in ionic strength. A collapse would only occur if the linker bends, 
and an observation of this phenomenon would support the solenoid model. 
The experiments by Yao et al. as well as more recent experiments by Butler 
and Thomas indeed reported a bending of the linkers but do not agree with 
experiments by Bednar et al. and by others that did not find any evidence 
for a collapse. 

In our work we shall follow the ideas put forward by Woodcock et al. [15] 
and Schiessel et al. [THJ but follow a rigorous mathematical approach. 

2 The two-angle-model 

Following Woodcock [TH] et al. and Schiessel et al. [IE] we consider four 
consecutive nucleosomes (cf. Fig. |TJ): N ,Ni,N 2 and N 3 e 1R 3 within the 
chain. N 3 is a function of N , ..,N 2 by fulfilling the following conditions: 
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z) <((N Q -N 1 ),(N 2 -N 1 )) = a; 

ii) ||iV 2 -JVi|| =62,11^0-^11 = fci, ||JV 3 — iV 2 || =63, with 61, ...,6 3 = 6; 

Hi) P := {r e R 3 | 3 A, ft e R, such that r = N x + X(N - JVk) + fi(N 2 - iVi)} 

P' := {r 6 I 3 I 3A', 11' e E, such that r = N± + X'(N 2 - N{) + /i'(N 3 - Nj} 

<(P,P')=f3. 

By straightforward considerations this leads to the following expression 
for iV 3 : 

(AT 2 - JVi) x (N - JVj) . N,-N 2 
' [|(iV 2 - JVO x (iVb - JV X )|| ' ' pVi-JVair 

where ||AT 2 — || = 6 2 arid Dt^, is the orthogonal rotational transformation 
matrix defined by the axis v e R 3 and the rotation angle (p e [0, 2n] (with 
respect to the right-hand- rule). (Note that the chromatin fibers described 
by these formula do not show a tangential distance between the ingoing and 
the outcoming DNA linkers). 

So the geometrical structure of the necklace is determined entirely by 
a, (3 and b. But this model only describes the linker geometry and does not 
account for any forms of nucleosome-nucleosome interaction. Furthermore it 
assumes straight linkers. It is still not completely clear whether the linkers 
of the 30nm fiber at high salt concentration are straight. But the bending of 
the linkers would cost about lOkT and straight linkers are very likely, cf. jjj. 

For every ideal chromatin fiber with a certain set of values (a, (3, b) 
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it is possible to construct a spiral with radius R and a gradient m so that 
all the nucleosomes are located on this spiral (cf. [19|). The nucleosomes 
are placed along the spiral in such a way that successive nucleosomes have a 
fixed (Euclidean) distance b from one another (in fact there are many such 
spirals, but the interesting spiral is the one with the largest gradient m). The 
parametrization of the spiral is given by 



^R-cos(^) { 



l{t) ■- 



t 



t e 



J 



where R is the radius and m — - is the gradient of this master solenoid 
(which follows from r(0) = (0, a, 1) T ). 

By straightforward calculations one can show the following relation: 



ir 2R 2 {l- cos ( ) + d 2 



R 

Furthermore, cos(7r — a) — ~^prp with r = -Ri — Rq, r = R 2 — R\ and 
r 2 = Rq — R-i leads to: 



cos 7r — a 



2R 2 cos (^) (1 - cos (^)) + d 2 



(2) 



2R? (1 -cos(^)) +d? 
Finally, j3 can be calculated by evaluating cos(/?) = 

angle between two successive planes (cf. iii): 

, d 2 cos + R 2 sin 2 (^) 

cos (0) = ^ 9/ „ KrJ 

y ! d 2 + R 2 sin 2 (^) 

So equations [TJ |21 and El relate R, a and d to a, (3 and b and thus the 
global fiber geometry to these local variables. These equations are important 
for further calculations. 



<roXn|r 2 xr > ^ 
||roXT-i||-||9~2 xj"o|| ' 



(3) 
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The inverse transformation, which relates the local fiber properties to 
the global spiral geometry (i.e. (a, (3, b) (R, m = -, d)) is very useful, too, 
and given by (cf. [IH]) : 

fe-oos(f) 
2(1 -sin 2 (f)cos* (§)) 



KJ \ 2 / 



m _l_ 2sin(f)^l-sin 2 (f)cos^(f) 
m a cot (f) arccos (2 sin 2 (f ) cos 2 (f ) - l) 



j= , 6sin(f) , + ■ w 

'esc 2 (f ) - cos 2 (f ) 2 V csc2 (f ) - 1 



These three equations relate the local fiber geometry to the global properties 
of the associated spiral. 



3 Phase Diagram 

Fig. |21 shows the chromatin phase diagram: For every point 

7„ e T n := {(ai)i=i... n , Oj G IR 3 | the series (a,) fulfills i)-iii) with a, (3 e [0, 7r]} 

there is a certain fiber structure (i.e. T n is the phase space of all ideal fibers 
of length n). These structures will be discussed in the following two sections. 
They can be divided into planar and three-dimensional structures. 
In fact the true phase diagram would be at least three-dimensional, be- 
cause the fiber's length N is an important parameter, too. It plays an 
important role from the numerical point, when considering different func- 
tions, which depend on the fiber's geometry as well as on its length such as 
nucleosome-density or excluded-volume-interactions. In these cases the fiber 
length should be as small as possible to lower the running time of the pro- 
gram but it should also be large enough to obtain the right physical insights. 
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This optimization principle is not always easy to fulfill - especially in the 
latter case of the hard-core excluded- volume-interactions. 
Another parameter which is not shown in Fig. |2]is the linker length b. But 
this can be easily ignored as the linker length is considered to be a constant 
for ideal fibers, and the mean value of it is fixed. 

The classical solenoid model of Finch and Klug (cf. is found in the large 
a and small (3 area in Fig. El (in their case the linkers were bent). Various 
other structures with a=30 and different values of (3 were found by Wood- 
cock et al. (cf. [T5]). These are located on the left side of Fig. |2j So all 
these possible 30nm-structures can be found in the phase diagram. 

If either one of the angles a or (3 equals or tt the resulting structure 
will be planar. Although it seems that the planar structures are not of much 
interest from the physical point of view it turned out that they play a very 
important role for the excluded volume structure of the phase diagram, and 
therefore they will be discussed here in detail. 

In the case (3 = tt and a arbitrary the fiber forms 2D zig-zag-like structures, 
as shown at the top of Fig. |21 The length of a fiber consisting of N monomers 
is given by L = d ■ N = b ■ sin (|) N (cf. Eq. EJ) and the diameter is given by 
D = b cos (^) (cf. Eq. 0}. The length of the fiber increases with increasing 
a. In the case of a = tt with an arbitrary value of (3 the fiber is a straight 
line and for a = and (3 e [0,tt] it corresponds to linkers that go back and 
forth between two positions (cf. Fig. EJ). 

Consider now the important case of the line (3 — at the bottom of Fig. |21 
If a ~ tt, the fiber forms a circle (with radius R ~ as follows from Eq. 
EJ). Its radius converges towards infinity for a — ► 7r. For certain values of a 
the fiber forms a regular polygon. For instance the value a = ~ corresponds 
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to the square and a = 60 is the regular triangle (cf. Fig. EJ). To characterize 
these regular polygons one needs two variables: At first the number of the 
tips i and secondly the order n of the polygon which gives the number of 
loops the fiber needs to arrive at the starting point again. 
These special values of a = af are given by 



The order n of a™ is a measure of its influence on the forbidden surfaces 
of the excluded volume structure of the phase diagram. Therefore all different 
values of a™ have to be sorted by n first and secondly by i (if ri\ > n 2 , then 
a™ 1 is "more important" than a™ 2 and if i\ > %2 then cu™ is "more important" 
than a") - this is what makes (A) necessary: For example a\ = a$ = ctg 
but the order of these three a™ is always n = 1 and therefore they are all in 
the same equivalence class (the special case i = 2n leads to one-dimensional 
structures with ct\ = - this value of a™ of highest order in n and i plays an 
important role for the forbidden area in the phase diagram and is therefore 
mentioned here, too). 

All af are irrational numbers (cf. Eq. Ej). For numerical reasons it is therefore 
useful to give all angles in degrees, because these are always rational numbers. 

Furthermore, (A) <^ $ a, j3,b E N with n = a ■ b and i = (3 ■ b, which 
means that n and i are coprime (it) as a™ depends only on the relation j. 

" =^>" Let n be an arbitrary element of N. Now assume that i G N, i > 2n and 
n and i have a common divisor b i.e. 3 a, (3 and b e N such that i = a-b and 




$ n', j'eN with n' < n and a" = a"' (A). 
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n = 8-b =>• 2 = = £ = i furthermore i > 2n indicates | > 2f =>• ^ < 

o 

I which means 3 n', i' e N with > 2n' and y = j, namely n' — (3 and = 
a. This means -i-fc implies -<A which is equivalent to A ^ 
"-<=" To prove the other implication above it suffices to show that if n and 
i are coprime, there are no numbers m, j with m < n such that nj = mi. 
Dividing this equation by gcd(n, m) one gets a new equation of the same 
structure with m! and n' instead of m and n, fulfilling gcd(m', n') = 1. 
Moreover, gcd(n, i) = 1 leads to gcd(n',i) = 1. Thus gcd(m'i,n') — 1 =>- 
j = 1 which contradicts to j > 2n. 

Now consider the case that n and n' are coprime with n' > n, i > 2n and 
%' > 2n' such that j, — £ j| (which means that they are not equivalent to 
any order- 1- value of a"). Furthermore, assume that — = | =>- «' = (^) • i, 
it is clear that (^) ^ N and ieN, therefore 3 a G N so that i = a-n which 
contradicts (-) ^ N. This means that if two orders ri\ and n 2 are coprime 
with 7ii,ri2 > 1 , there are never equal values of a™ 1 and a™ 2 for all possible 
« > 2n, i' > 2n' and «, i' e N. For example, n and n + 1 are always coprime 
numbers and therefore have never common a™. 

So for a given order n the possible values of af depend on the prime 
factor dismantling of n and therefore are very irregular. If n is a prime 
number, all i > 2n with % £ n ■ N are allowed: for n = 2 alH > 4 with 
i £ 2 ■ N are allowed and for n = 3 all i with i > 2n, i ^ 3 • N are allowed 
(and for n — 1 alH > 2 are allowed). So between two values of a* there is 
one of of and two of af . 

The distances A™ = a™ — a™ between two consecutive i 2 > %\ are given 
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below for the first three and additional for all prime orders: 



n = 1 : a- = tt -, (i > 2n) =>- A, = — — 

9 47T . , . 9 87T 

n = 2 : a, 2 = tt - — , (i > 2n and i 2 ■ N) =>- A 2 - 



i ' v " " 7 i(i + 2) 

- 6ff ^o„„„,^ 9 m^a3_ mod(*,3)-67r 



n = 3 : af = n , (i>2n and i 4 3 • N) A, 

2 v ~ 7 ? ' i(« + mod 0,3)) 

n ■ 2n . . . A „ m • n • 27T 
n prime : a, = 7r , u > 2n and iff n • N) =>- A • = — . 

i i(i + m) 



This shows furthermore that A" — > 0, because the counter is always ~ n 
and the leading term of the denominator is at least ~ n 2 (as % > 2n). 

The most interesting structures of the chromatin phase diagram are the 
following two cases: solenoid-like structures and fibers with crossed-linkers. 
An extensive discussion of these two structures can be found in ^H] and here 
a short overview will be given. 

For small (3 <C 1 and a ~ tt the chromatin fibers resemble solenoids 
where the linkers themselves follow closely a helical path. This leads to 
the condition \jjf) <C 1. The lowest order approximation in (^) indicates 
b » dv/TT^ 2 (cf. Eq. CD, a « vr - 7 ^ 5 y (cf. Eq. ED and (3 » ^ 
(cf. Eq. ED- These approximations lead to the geometrical properties of the 
solenoid-like fibers: The radius of the fiber is r, is the length of a fiber of 
N monomers, A = N/In is the line density and g = X/irr 2 is the density of 
the fiber: 



b (tt - a) bN (3 

(3 2 + (<k- of W "^Tp? 

y//3 2 + (tt - cQ 2 = (/j 2 + (tt - a) 2 ) 5 / 2 

bp Q ~ n b 3 [3 (n - a) 2 
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The vertical distance S between two loops plays an important role in the 
following sections. It can be obtained by 5 ~ I 2* , which leads to 

tv — a. 

(-*-) b(3 

S « - u ~ aJ = . (8) 
One can calculate the exact value of 5 by: 

27rr/a 

5= | m( 7 (o)-i(y,(o,y tf (*)) t i* 

467rsin(f)sin(f) 
arccos (2 cos 2 (|) sin 2 (|) — l) a/3 + cos(a) + (cos(a) — 1) cos(/3) 

However, the approximation (jHJ) is much more useful. 

For p <C (7T — a) one finds a very dense spiral (5 <C r) and in the opposite 
limit (3 ^> (tt — a) the solenoid has a very open structure >r). 

Next we turn to the crossed-linker fibers. Now consider the case /3 <C 1 
and a 7r: Above the regular polygons were discussed. For a non-vanishing 
/? these regular polymers open up in an accordion-like manner. This leads to 
a three-dimensional fiber with crossed linkers (cf. Fig. |2J. Using and 
Eq. □ Eq. EJand Eq. El one gets (cf. US]): 

d 2 = ^(3 2 {AR 2 -b 2 ) 



and thus 



(3 2 ^r 77 ~~ a 



1 cot 



2cos(a/2) V 4 2 

N j3 b 7r — a 
1 n = cot(^— ) 

4 16 



Q 



a fib' ir (3 a b 3 

Now 5 follows again from S ~ I 2^ : 
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Of course, not all points 7 n G T n have the same probability to occur in 
real chromatin fibers. There are even some forbidden areas in the phase space 
due to excluded volume restrictions. These forbidden areas are the subject 
of this section. In this numerical and analytical study hard-core potentials 
were used to model the excluded volume interactions of the fibers. The border 
which separates the non-excluded volume structures from those which fulfill 
the excluded volume condition is denoted by ( (which is a function of a and 
is plotted in Fig. El too). 

n m i n = 3 (-3 — l) = 3 f 2 ^) is given by the chromatin geometry (linker 
length, octamer diameter and height, see below). As shown in Fig. |2]the in- 
teresting part of the phase diagram for the excluded-volume phase-transition 
is the lower one with (5 e [0, 30]. This cut-out of the phase diagram is shown 
in Fig. El 

One can roughly explain Fig. |3] by dividing the excluded volume inter- 
actions between the nucleosomes in short-range and long-range interactions: 
Between two consecutive nucleosomes there is never an excluded volume in- 
teraction, short-range interactions occur between nucleosome k and nucleo- 
some k + 2, and long-range interactions are those between nucleosomes with 
a larger distance than 2 (Ak > 2). These approximations can be found in 
HU, too. 

Another explanation of the structure of the excluded-volume phase- 
transition, which even allows to derive the explicit behavior of ( and does 
not need to divide into short-range and long-range interactions, is presented 
below. 

Every peak of the ^-function corresponds to a regular polygon: The large 
peaks correspond to polygons of order 1 and the smaller ones to order 2 and 
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3 polygons. Between two order-l-peaks there is one order-2-peak and two 
order-3-peaks. The classification of these £-peaks is shown in Fig. |3 the (- 
function has local maxima at every a™. The planar structures which belong 
to (af,3 = 0) need a large rise of 3 to arrive in the area of the excluded 
volume structures, because at 3 = nucleosome k and k + n ' 2n n are located 

' " 7T— af 

exactly at the same position. 

At first, consider only the special values a™: At first order increasing 3 from 
3 — to (3 — (3 leads to a vertical movement A™(/3) of the relevant micro- 
somes k and k+ n '^\ along the chromatin axis . If A™(/3) = d =>- C( a ?) = $■ 
For large a, Af(3) can be calculated by Eq. Eland one finds: 

Am = v iJ (io) 

with 3 > this leads to 

' (Ar(7r - o) 2 



and furthermore [3f(d) = ((af) implies 



'IF 



n\2 



cw)= rf2 ; (T "^'' • < 12 > 

Fig. 13 shows the numerically calculated zeta-function and the theoretical 
predictions for the maxima. As Eq. ITU1 shows. A"(/3) ~ n. This means that 
planar structures at an a™ of higher order need a smaller rise of 3 to fulfill 
the excluded volume conditions and therefore ( is decreasing with increasing 
a. In fact this is the reason why n is called the order of the peaks of the 
C-function. It can be understood easily if one considers the fact that a higher 
order means more nucleosomes are located between two overlapping ones. 
And thus a rise of 3 has a stronger effect than at lower orders. 
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As A™(/3) — > for a™ — > 7r the maxima of the ^-function converge to- 
wards for a™ — > 7r. There are infinitely many a" for n — > oo and the 
distance A™ converges to zero, so the ^-function has infinitely many local 
maxima and minima for a" — > it. It is clear that ((rr) = as the fiber forms 
a fully stretched fiber (a circle of radius r = oo). This explains the forbidden 
strip at the left side of the figures as the maximum of highest order (1,2) 
with a corresponding one-dimensional structure. 

So far the ^-function is only known at the positions of the maxima af. 
Now consider values of a, which are close to an a™, say a' = «f±Aa ("close" 
means such Act, which lead to a shift Aa < 2r). This leads to a slight shift 
Aa of those nucleosomes which are located at the same places (namely k and 
k+ (^r^))- At P — this shift Aa is orthogonal to the fiber's axis. This time 
A" still denotes the distance between the nucleosomes k and k+ (^r^) along 
the fiber's axis, but now they are not located at the same spots but slightly 
shifted. Therefore their distance A is no longer equal to their distance A" 
along the axis, when increasing (3. 

In this case increasing (3 still leads to a movement along the vertical axis of 
the fiber, but now the distance A of nucleosome k and k + (^zl) increases 
like A 2 m Aa 2 + (A™) 2 . The fiber fulfills excluded volume, when A = 2r, 
which means 

A^ = V4r 2 - Aa 2 < 2r. 

So the critical value of A™, which has to be achieved to fulfill excluded vol- 
ume, decreases with increasing Aa. Therefore the C-fimction has a local 
maximum at a": £(a:" ± Aa) < C( a i)- 
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To calculate Aai jn (Aa) imagine a planar structure of j nucleosomes with 
an entry-exit-angle a' = a™ ± Act of two consecutive octamers. The locations 
of these nucleosomes are denoted by p , pi, p 2 , Pj-\ G M 3 . Without loss 
of generality one can assume 



Po 







v / 



Pi 



b cosffl 



b sin(f) 



V 



and 



/ 



Wk>2 : p k = p k ^i + K ■ - p k - 2 ) 





( ( >\ 


sin(7r — a') 






cos(7r — a ) 




with ft := Ul_ al = 


— sin(7r — a') 


cos(7r — a') 







K 





V 



the rotational matrix along the z-axis. This leads to 

k-i 



Pk 



J2^ m -Pt V k > 0. 



m=0 



Now Adi t n is given by Aai in = ||pfc=i||, which leads to 



Aa,-, 



m=l 



(13) 



(were TZ and p\ depend on i, n and Aa) Now one can use Eq.^Jto calculate 
the ^-function around the maximal values a": 



v/C(a') 2 + (t - "0 s 



= \/4r 2 - Aa 2 



C(a = a? ± Aa) 



(4r 2 -|| Y. n ' n ' 1 -Pi || 2 )(vr-a') 2 

m=l 

\ " '' m=l 



(14) 



Similar to the case of large a one can find equations for the (^-function of 
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small a: Eq. El leads to 



with a' = a" ± Aa. Now again ((a') follows from /3(A? = * Ar 2 — Aaf n ) 



C(a'): 



C(a' = a? ± Aa) 



cot(f) 



n 



4f 2 



m— 1 . || 2 



(15) 



m=l 



|Aa| < c fulfilling Aa(c) = 2r gives the interval of the allowed a- values: 
a' G [af — c, a™ + c] for a certain peak (n,i). The predictions of Eqs. Upland 
ITBI are shown in Fig. El together with the simulation results. 



4 Excluded volume restrictions of the DNA 
linkers 

As mentioned before the nucleosome-nucleosome excluded volume interac- 
tions are not the only ones within the chromatin fiber. The DNA linkers 
have a diameter of about 2nm and therefore excluded volume restrictions, 
too. This is in particular very important for all crossed-linker structures. But 
since the DNA linkers have a very strong (although screened) Coulomb re- 
pulsion their excluded volume interactions can be revealed by looking at the 
Coulomb energy between the linkers. Other potentials like the nucleosome- 
nucleosome interaction or the interaction between DNA linkers and nucleo- 
somes will be neglected here. 

One can use the Debye-Hiickel theory to model the Coulomb repulsion of 
the DNA linkers, but since the screening of this interaction starts at the ra- 
dius of the DNA strand and due to the condensation of ions along the DNA 
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linkers, one has to calculate a correction of the screened potential by fitting 
the tail of the Debye-Hiickel potential for an infinitely long cylinder to the 
Gouy-Chapman potential in the far zone. This calculation can be found in 
fK)\ l2*T] and leads to a corrected linear charge density v e s which is also given 
in the table below for different levels of monovalent salt concentration and 
can be found in [22] for instance. During the last two to three decades DNA 
models based on such potentials have been developed and applied widely, 
and their predictions are usually very good [23 121123 123 123 123 HH • 



c s [10- 2 M] 


1 


2 


3 


4 


5 


k [nm _1 ] 


0.330 


0.467 


0.572 


0.660 


0.738 


u cS [e/nm] 


-2.43 


-2.96 


-3.39 


-3.91 


-4.15 



So one can calculate the Coulomb repulsion between the DNA segments i 
and j by evaulating 

= — / dxidxj (16) 

c J J r i,j 

where c is the total dielectric constant of water. These two integrals were 
numerically calculated and the results for the Coulomb energy of a single 
chromatin linker within a fiber is shown in Fig. 

One can see that the Coulomb repulsion of the linkers is very high within 
the gaps of the excluded volume borderline. The repulsion also diverges for 
the crossed linker fibers when f3 becomes too small. 

5 Discussion and Conclusion 

The following effects have to be considered to explain the differences between 
the calculated theoretical predictions and the simulational results. 
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The effective entry-exit angle a e is the projection of a onto a plane which 
is orthogonal to the axis of the master solenoid. It decreases with increasing 
(3. This is a consequence of the fact that the length of the fiber increases 
with increasing (3. a e is important for the calculation of the number of linkers 
which form a closed loop: Ni = 27T ( 2w gives the number of breaks of a e 
which has to be done to get a full loop. If is an integer, this is equal to 
the number of linkers, which corresponds to a full loop. If it is not an integer 
number, then the fractional part gives the fraction of 7r — a e which is missing 
for a full loop, but the entry-exit angle is fixed here and it is assumed that 
Ni = ). In the calculations above a e = a =const. was assumed, because 
only small (3 were considered. a e will be calculated as a function of a and (3 
below. It converges towards for (3 — *> 7r for different values of a and f3. 
To calculate a e consider three consecutive nucleosome locations within a 
given fiber: n , n\ and n<i el 3 . Without loss of generality one can assume: 
ra o = Po, Tit = pi + po and n 2 = P2 + Pi + Po were the pi are the following 
linker vectors: 



Po 







v / 



^Vb 2 - d 2 ^ 



Pi 



\ 





d 



and p 2 



J 



( \ 



This means that the z-axis of the coordinate system is the axis of the fiber 
(d is given by Eq. Inland b is the linker length). Now ||p 2 || = b leads to 



b 2 -d 2 



and %2 can be calculated by using 



cos a 



x 2 



_ (-Pi I P2) 
b 2 

b 2 cos(a) + d 2 
- d 2 
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with a = —, = = ' — == this leads to 

v / csc 2( f )_ cos 2(|) 



, (b 2 cos(a) + d 2 ) 2 d /,„ „ f /3\ , 2 . 
y 2 = \jb 2 -d 2 - - b 2_ d 2 — = \i h cos I o ) sm 2 (a) 



and now a e follows from the projection onto the x-?/-plane: 



sin 



a. 



2^/x 2 + y 2 

sin (^ )= >(f)^;f . (17) 

1 2 J V 2(l + cos(a)) V ; 

As Eq. IT7I shows. a e decreases with increasing (3. As a consequence of this 
Ni decreases, too: If a e — > 0, then Ni(a,j3) = 7r _ Q 2 | Q , ^ — > 2. So all fibers 
with high values of beta need only approx. 2 nucleosomes for a complete 
loop. This is confirmed by Fig. |21 As a e was assumed to be constant, Ni 
was also constant in the calculations above. This is a suitable approximation 
for small f3. For larger f3 the assumed values of Ni were a bit too large and 
therefore the calculated values of ( were a bit too small - but this effect is 
small compared to the other estimations which were made above. The error 
of (, due to the assumption that Ni is constant, increases with increasing f3. 



So to increase accuracy in the calculations above one should replace the 

term -^7 by 2 7 , m and furthermore one should use Eq. [HI to derive exact 

equations in order to replace Eqs. |H1 and But this can only be done by 
using numerical methods and not analytical ones. Moreover, Fig. |H1 shows 
that the assumed approximations lead to useful equations, which describe ( 
in a suitable way. 
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Figure 2: The chromatin phase diagram with some chromatin structures for 
different values of a and (3. The solenoid and crossed-linker structures are 
most important. The dotted line is the function ((a) which is explained in 
the text. 
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chromatin phase diagram - excluded volume "phase transition" 
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Figure 3: Fine-structure of the excluded- volume "phase transition". The 
chromatin fibers in the green area fulfill the excluded volume conditions, 
those in the yellow area do not. The borderline is the function C{ a )- 
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excluded volume phase transition and its theoretical prediction 




Figure 4: Calculated theoretical prediction (red) and the simulation result 
(blue) for ( and small a. 
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excluded volume phase transition and its theoretical prediction 
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Figure 5: Calculated theoretical prediction (red) and the simulation result 
(blue) for ( and larger a. 
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The ^-function : classification of its maxima 
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Figure 6: Classification of the excluded volume regions. 
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Coulomb repulsion of the DNA linkers 
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Figure 7: Coulomb repulsion of the DNA linkers. Shown are also two cuts 
in the inset. 
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